library(mosaic)
Loading required package: dplyr
Registered S3 method overwritten by 'dplyr':
  method           from
  print.rowwise_df     

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

Loading required package: lattice
Loading required package: ggformula
Loading required package: ggplot2
Loading required package: ggstance

Attaching package: ‘ggstance’

The following objects are masked from ‘package:ggplot2’:

    geom_errorbarh, GeomErrorbarh


New to ggformula?  Try the tutorials: 
    learnr::run_tutorial("introduction", package = "ggformula")
    learnr::run_tutorial("refining", package = "ggformula")
Loading required package: mosaicData
Loading required package: Matrix
Registered S3 methods overwritten by 'htmltools':
  method               from         
  print.html           tools:rstudio
  print.shiny.tag      tools:rstudio
  print.shiny.tag.list tools:rstudio
Registered S3 method overwritten by 'htmlwidgets':
  method           from         
  print.htmlwidget tools:rstudio
Registered S3 method overwritten by 'mosaic':
  method                           from   
  fortify.SpatialPolygonsDataFrame ggplot2

The 'mosaic' package masks several functions from core packages in order to add 
additional features.  The original behavior of these functions should not be affected by this.

Note: If you use the Matrix package, be sure to load it BEFORE loading mosaic.

Attaching package: ‘mosaic’

The following object is masked from ‘package:Matrix’:

    mean

The following object is masked from ‘package:ggplot2’:

    stat

The following objects are masked from ‘package:dplyr’:

    count, do, tally

The following objects are masked from ‘package:stats’:

    binom.test, cor, cor.test, cov, fivenum, IQR, median, prop.test, quantile,
    sd, t.test, var

The following objects are masked from ‘package:base’:

    max, mean, min, prod, range, sample, sum
library(tidyverse)
Registered S3 methods overwritten by 'dbplyr':
  method         from
  print.tbl_lazy     
  print.tbl_sql      
── Attaching packages ───────────────────────────────────────────────── tidyverse 1.3.0 ──
✓ tibble  2.1.3     ✓ purrr   0.3.3
✓ tidyr   1.0.0     ✓ stringr 1.4.0
✓ readr   1.3.1     ✓ forcats 0.4.0
── Conflicts ──────────────────────────────────────────────────── tidyverse_conflicts() ──
x mosaic::count()            masks dplyr::count()
x purrr::cross()             masks mosaic::cross()
x mosaic::do()               masks dplyr::do()
x tidyr::expand()            masks Matrix::expand()
x dplyr::filter()            masks stats::filter()
x ggstance::geom_errorbarh() masks ggplot2::geom_errorbarh()
x dplyr::lag()               masks stats::lag()
x tidyr::pack()              masks Matrix::pack()
x mosaic::stat()             masks ggplot2::stat()
x mosaic::tally()            masks dplyr::tally()
x tidyr::unpack()            masks Matrix::unpack()
library(lubridate)

Attaching package: ‘lubridate’

The following object is masked from ‘package:base’:

    date
library(DataComputing)
library(rvest)
Loading required package: xml2

Attaching package: ‘rvest’

The following object is masked from ‘package:purrr’:

    pluck

The following object is masked from ‘package:readr’:

    guess_encoding

Guiding Question (i.e., Research Question):

How do different aspects of determining how developed a country is (such as birth rate, population, death rate, etc.) manifest themselves in the spread of COVID-19. Essentially, how does COVID-19 show different levels of impact on different countries?

Data Access

Reading in the Data:

Data Source 1: COVID

COVID <- read.csv(file = "total-covid-cases-deaths-per-million.csv")
COVID
COVID %>%
  nrow()
[1] 9487
COVID %>%
  names()
  [1] "total.covid.cases.deaths.per.million" "X"                                   
  [3] "X.1"                                  "X.2"                                 
  [5] "X.3"                                  "X.4"                                 
  [7] "X.5"                                  "X.6"                                 
  [9] "X.7"                                  "X.8"                                 
 [11] "X.9"                                  "X.10"                                
 [13] "X.11"                                 "X.12"                                
 [15] "X.13"                                 "X.14"                                
 [17] "X.15"                                 "X.16"                                
 [19] "X.17"                                 "X.18"                                
 [21] "X.19"                                 "X.20"                                
 [23] "X.21"                                 "X.22"                                
 [25] "X.23"                                 "X.24"                                
 [27] "X.25"                                 "X.26"                                
 [29] "X.27"                                 "X.28"                                
 [31] "X.29"                                 "X.30"                                
 [33] "X.31"                                 "X.32"                                
 [35] "X.33"                                 "X.34"                                
 [37] "X.35"                                 "X.36"                                
 [39] "X.37"                                 "X.38"                                
 [41] "X.39"                                 "X.40"                                
 [43] "X.41"                                 "X.42"                                
 [45] "X.43"                                 "X.44"                                
 [47] "X.45"                                 "X.46"                                
 [49] "X.47"                                 "X.48"                                
 [51] "X.49"                                 "X.50"                                
 [53] "X.51"                                 "X.52"                                
 [55] "X.53"                                 "X.54"                                
 [57] "X.55"                                 "X.56"                                
 [59] "X.57"                                 "X.58"                                
 [61] "X.59"                                 "X.60"                                
 [63] "X.61"                                 "X.62"                                
 [65] "X.63"                                 "X.64"                                
 [67] "X.65"                                 "X.66"                                
 [69] "X.67"                                 "X.68"                                
 [71] "X.69"                                 "X.70"                                
 [73] "X.71"                                 "X.72"                                
 [75] "X.73"                                 "X.74"                                
 [77] "X.75"                                 "X.76"                                
 [79] "X.77"                                 "X.78"                                
 [81] "X.79"                                 "X.80"                                
 [83] "X.81"                                 "X.82"                                
 [85] "X.83"                                 "X.84"                                
 [87] "X.85"                                 "X.86"                                
 [89] "X.87"                                 "X.88"                                
 [91] "X.89"                                 "X.90"                                
 [93] "X.91"                                 "X.92"                                
 [95] "X.93"                                 "X.94"                                
 [97] "X.95"                                 "X.96"                                
 [99] "X.97"                                 "X.98"                                
[101] "X.99"                                 "X.100"                               
[103] "X.101"                                "X.102"                               
[105] "X.103"                                "X.104"                               
[107] "X.105"                                "X.106"                               
[109] "X.107"                                "X.108"                               
[111] "X.109"                                "X.110"                               
[113] "X.111"                                "X.112"                               
[115] "X.113"                                "X.114"                               
[117] "X.115"                                "X.116"                               
[119] "X.117"                                "X.118"                               
[121] "X.119"                                "X.120"                               
[123] "X.121"                                "X.122"                               
[125] "X.123"                                "X.124"                               
[127] "X.125"                                "X.126"                               
[129] "X.127"                                "X.128"                               
[131] "X.129"                                "X.130"                               
[133] "X.131"                                "X.132"                               
[135] "X.133"                                "X.134"                               
[137] "X.135"                                "X.136"                               
[139] "X.137"                                "X.138"                               
[141] "X.139"                                "X.140"                               
[143] "X.141"                                "X.142"                               
[145] "X.143"                                "X.144"                               
[147] "X.145"                                "X.146"                               
[149] "X.147"                                "X.148"                               
[151] "X.149"                                "X.150"                               
[153] "X.151"                                "X.152"                               
[155] "X.153"                                "X.154"                               
[157] "X.155"                                "X.156"                               
[159] "X.157"                                "X.158"                               
[161] "X.159"                                "X.160"                               
[163] "X.161"                                "X.162"                               
[165] "X.163"                                "X.164"                               
[167] "X.165"                                "X.166"                               
[169] "X.167"                                "X.168"                               
[171] "X.169"                                "X.170"                               
[173] "X.171"                                "X.172"                               
[175] "X.173"                                "X.174"                               
[177] "X.175"                                "X.176"                               
[179] "X.177"                                "X.178"                               
[181] "X.179"                                "X.180"                               
[183] "X.181"                                "X.182"                               
[185] "X.183"                                "X.184"                               
[187] "X.185"                                "X.186"                               
[189] "X.187"                                "X.188"                               
[191] "X.189"                                "X.190"                               
[193] "X.191"                                "X.192"                               
[195] "X.193"                                "X.194"                               
[197] "X.195"                                "X.196"                               
[199] "X.197"                                "X.198"                               
[201] "X.199"                                "X.200"                               
[203] "X.201"                                "X.202"                               
[205] "X.203"                                "X.204"                               
[207] "X.205"                                "X.206"                               
[209] "X.207"                                "X.208"                               
[211] "X.209"                                "X.210"                               
[213] "X.211"                                "X.212"                               
[215] "X.213"                                "X.214"                               
[217] "X.215"                                "X.216"                               
[219] "X.217"                                "X.218"                               
[221] "X.219"                                "X.220"                               
[223] "X.221"                                "X.222"                               
[225] "X.223"                                "X.224"                               
[227] "X.225"                                "X.226"                               
[229] "X.227"                                "X.228"                               
[231] "X.229"                                "X.230"                               
[233] "X.231"                                "X.232"                               
[235] "X.233"                                "X.234"                               
[237] "X.235"                                "X.236"                               
[239] "X.237"                                "X.238"                               
[241] "X.239"                                "X.240"                               
[243] "X.241"                                "X.242"                               
[245] "X.243"                                "X.244"                               
[247] "X.245"                                "X.246"                               
[249] "X.247"                                "X.248"                               
[251] "X.249"                                "X.250"                               
[253] "X.251"                                "X.252"                               
[255] "X.253"                                "X.254"                               
COVID %>%
  head()

Data Source 2: CountryData

CountryData
CountryData %>%
  nrow()
[1] 256
CountryData %>%
  names()
 [1] "country"           "area"              "pop"               "growth"           
 [5] "birth"             "death"             "migr"              "maternal"         
 [9] "infant"            "life"              "fert"              "health"           
[13] "HIVrate"           "HIVpeople"         "HIVdeath"          "obesity"          
[17] "underweight"       "educ"              "unemploymentYouth" "GDP"              
[21] "GDPgrowth"         "GDPcapita"         "saving"            "indProd"          
[25] "labor"             "unemployment"      "family"            "tax"              
[29] "budget"            "debt"              "inflation"         "discount"         
[33] "lending"           "narrow"            "broad"             "credit"           
[37] "shares"            "balance"           "exports"           "imports"          
[41] "gold"              "externalDebt"      "homeStock"         "abroadStock"      
[45] "elecProd"          "elecCons"          "elecExp"           "elecImp"          
[49] "elecCap"           "elecFossil"        "elecNuc"           "elecHydro"        
[53] "elecRenew"         "oilProd"           "oilExp"            "oilImp"           
[57] "oilRes"            "petroProd"         "petroCons"         "petroExp"         
[61] "petroImp"          "gasProd"           "gasCons"           "gasExp"           
[65] "gasImp"            "gasRes"            "mainlines"         "cell"             
[69] "netHosts"          "netUsers"          "airports"          "railways"         
[73] "roadways"          "waterways"         "marine"            "military"         
CountryData %>%
  head()

Data Source 1: Continents

Continents <- read.csv(file = "countries and continents.csv")
Continents
Continents %>%
  nrow()
[1] 251
Continents %>%
  names()
 [1] "name"                             "official_name_en"                
 [3] "official_name_fr"                 "ISO3166.1.Alpha.2"               
 [5] "ISO3166.1.Alpha.3"                "M49"                             
 [7] "ITU"                              "MARC"                            
 [9] "WMO"                              "DS"                              
[11] "Dial"                             "FIFA"                            
[13] "FIPS"                             "GAUL"                            
[15] "IOC"                              "ISO4217.currency_alphabetic_code"
[17] "ISO4217.currency_country_name"    "ISO4217.currency_minor_unit"     
[19] "ISO4217.currency_name"            "ISO4217.currency_numeric_code"   
[21] "is_independent"                   "Capital"                         
[23] "Continent"                        "TLD"                             
[25] "Languages"                        "Geoname.ID"                      
[27] "EDGAR"                           
Continents %>%
  head()

Data Wrangling of COVID Dataset

COVID
TidyCOVID <- COVID %>%
  rename(c("country" = "total.covid.cases.deaths.per.million") ) %>%
  rename(c("Code" = "X") ) %>%
  rename(c("Date" = "X.1") ) %>%
  rename(c("DeathsPerMillion" = "X.2")) %>%
  rename(c("CasesPerMillion" = "X.3")) %>%
  filter(row_number() > 1) %>%
  subset(select = c(1,2,3,4,5)) %>%
  mutate( country = as.character(country) ) %>%
  mutate( Code = as.character(Code) ) %>%
  mutate(Date = mdy(Date)) %>%
  mutate(DeathsPerMillion = as.integer(DeathsPerMillion) - 1) %>%
  mutate(CasesPerMillion = as.integer(CasesPerMillion) - 1)
Error: All arguments must be named

Data Wrangling of Continents

Data Visualization

JOEY: I think this is a good one to start out with, but if there is some way that you could make this more clear (color code, etc.) that would be great. This is a good start to our data.

JOEY: I like this graph because it visualizes one of the factors that contributes to the development classification. Can we add in a color code or something that would represent one more factor, such as different countries or regions of the world. The more that we can do with this one, the better.

JOEY: I also liked this graph, but we need to remove the outlier that scews it and makes the scale hard to see a relationship. It is the one with an X-value of over 10,000

JOEY: I think one more graph here will be sufficient. Can we do one that compares the relationship between Cases per Million, Deaths per million, and population!!! Maybe something like a cloropleth graph (I think thats what I’m thinking of)?

LS0tCnRpdGxlOiAiRmluYWwgUHJvamVjdCBSb3VnaCBEcmFmdCIKb3V0cHV0OiBodG1sX25vdGVib29rCmF1dGhvcnM6ICJFdmVseW4gTXVycmF5IGFuZCBKb3NlcGggUGV2bmVyIgotLS0KCmBgYHtyfQpsaWJyYXJ5KG1vc2FpYykKbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkobHVicmlkYXRlKQpsaWJyYXJ5KERhdGFDb21wdXRpbmcpCmxpYnJhcnkocnZlc3QpCmBgYAoKIyMgR3VpZGluZyBRdWVzdGlvbiAoaS5lLiwgUmVzZWFyY2ggUXVlc3Rpb24pOgoKSG93IGRvIGRpZmZlcmVudCBhc3BlY3RzIG9mIGRldGVybWluaW5nIGhvdyBkZXZlbG9wZWQgYSBjb3VudHJ5IGlzIChzdWNoIGFzIGJpcnRoIHJhdGUsIHBvcHVsYXRpb24sIGRlYXRoIHJhdGUsIGV0Yy4pIG1hbmlmZXN0IHRoZW1zZWx2ZXMgaW4gdGhlIHNwcmVhZCBvZiBDT1ZJRC0xOS4gIEVzc2VudGlhbGx5LCBob3cgZG9lcyBDT1ZJRC0xOSBzaG93IGRpZmZlcmVudCBsZXZlbHMgb2YgaW1wYWN0IG9uIGRpZmZlcmVudCBjb3VudHJpZXM/CgoKIyMgRGF0YSBBY2Nlc3MKClJlYWRpbmcgaW4gdGhlIERhdGE6CgoKRGF0YSBTb3VyY2UgMTogQ09WSUQKYGBge3J9CkNPVklEIDwtIHJlYWQuY3N2KGZpbGUgPSAidG90YWwtY292aWQtY2FzZXMtZGVhdGhzLXBlci1taWxsaW9uLmNzdiIpCkNPVklECmBgYAoKYGBge3J9CkNPVklEICU+JQogIG5yb3coKQpgYGAKYGBge3J9CkNPVklEICU+JQogIG5hbWVzKCkKYGBgCmBgYHtyfQpDT1ZJRCAlPiUKICBoZWFkKCkKYGBgCgoKCkRhdGEgU291cmNlIDI6IENvdW50cnlEYXRhCmBgYHtyfQpDb3VudHJ5RGF0YQpgYGAKCmBgYHtyfQpDb3VudHJ5RGF0YSAlPiUKICBucm93KCkKYGBgCmBgYHtyfQpDb3VudHJ5RGF0YSAlPiUKICBuYW1lcygpCmBgYApgYGB7cn0KQ291bnRyeURhdGEgJT4lCiAgaGVhZCgpCmBgYAoKCgpEYXRhIFNvdXJjZSAxOiBDb250aW5lbnRzCmBgYHtyfQpDb250aW5lbnRzIDwtIHJlYWQuY3N2KGZpbGUgPSAiY291bnRyaWVzIGFuZCBjb250aW5lbnRzLmNzdiIpCkNvbnRpbmVudHMKYGBgCmBgYHtyfQpDb250aW5lbnRzICU+JQogIG5yb3coKQpgYGAKYGBge3J9CkNvbnRpbmVudHMgJT4lCiAgbmFtZXMoKQpgYGAKYGBge3J9CkNvbnRpbmVudHMgJT4lCiAgaGVhZCgpCmBgYAoKCgoKIyMjIERhdGEgV3JhbmdsaW5nIG9mIENPVklEIERhdGFzZXQKCmBgYHtyfQpDT1ZJRApgYGAKCgpgYGB7cn0KVGlkeUNPVklEIDwtIENPVklEICU+JQogIHJlbmFtZShjKCJjb3VudHJ5IiA9ICJ0b3RhbC5jb3ZpZC5jYXNlcy5kZWF0aHMucGVyLm1pbGxpb24iKSApICU+JQogIHJlbmFtZShjKCJDb2RlIiA9ICJYIikgKSAlPiUKICByZW5hbWUoYygiRGF0ZSIgPSAiWC4xIikgKSAlPiUKICByZW5hbWUoYygiRGVhdGhzUGVyTWlsbGlvbiIgPSAiWC4yIikpICU+JQogIHJlbmFtZShjKCJDYXNlc1Blck1pbGxpb24iID0gIlguMyIpKSAlPiUKICBmaWx0ZXIocm93X251bWJlcigpID4gMSkgJT4lCiAgc3Vic2V0KHNlbGVjdCA9IGMoMSwyLDMsNCw1KSkgJT4lCiAgbXV0YXRlKCBjb3VudHJ5ID0gYXMuY2hhcmFjdGVyKGNvdW50cnkpICkgJT4lCiAgbXV0YXRlKCBDb2RlID0gYXMuY2hhcmFjdGVyKENvZGUpICkgJT4lCiAgbXV0YXRlKERhdGUgPSBtZHkoRGF0ZSkpICU+JQogIG11dGF0ZShEZWF0aHNQZXJNaWxsaW9uID0gYXMuaW50ZWdlcihEZWF0aHNQZXJNaWxsaW9uKSAtIDEpICU+JQogIG11dGF0ZShDYXNlc1Blck1pbGxpb24gPSBhcy5pbnRlZ2VyKENhc2VzUGVyTWlsbGlvbikgLSAxKQoKVGlkeUNPVklECgpgYGAKCgpgYGB7cn0KVGlkeUNPVklECgpgYGAKCiMjIyBEYXRhIFdyYW5nbGluZyBvZiBDb250aW5lbnRzCmBgYHtyfQpDb250aW5lbnRzCgpgYGAKCmBgYHtyfQpUaWR5Q29udCA8LQogIENvbnRpbmVudHMgJT4lCiAgZmlsdGVyKHJvd19udW1iZXIoKSA+IDIpICU+JQogIHN1YnNldChzZWxlY3QgPSBjKDEsMjMpKSAlPiUKICByZW5hbWUoYygiY291bnRyeSIgPSAibmFtZSIpICkgJT4lCiAgbXV0YXRlKCBjb3VudHJ5ID0gYXMuY2hhcmFjdGVyKGNvdW50cnkpICkgJT4lCiAgbXV0YXRlKCBDb250aW5lbnQgPSBhcy5jaGFyYWN0ZXIoQ29udGluZW50KSApCgpUaWR5Q29udAoKCmBgYAoKCgojIyBEYXRhIFZpc3VhbGl6YXRpb24KYGBge3J9CkRhdGFNaXggPC0gbGVmdF9qb2luKFRpZHlDT1ZJRCwgQ291bnRyeURhdGEsIGJ5ID0gImNvdW50cnkiKQpgYGAKCkpPRVk6IEkgdGhpbmsgdGhpcyBpcyBhIGdvb2Qgb25lIHRvIHN0YXJ0IG91dCB3aXRoLCBidXQgaWYgdGhlcmUgaXMgc29tZSB3YXkgdGhhdCB5b3UgY291bGQgbWFrZSB0aGlzIG1vcmUgY2xlYXIgKGNvbG9yIGNvZGUsIGV0Yy4pIHRoYXQgd291bGQgYmUgZ3JlYXQuICBUaGlzIGlzIGEgZ29vZCBzdGFydCB0byBvdXIgZGF0YS4KYGBge3J9CgpnZ3Bsb3QoZGF0YSA9IERhdGFNaXgsYWVzKHggPSBEYXRlLCB5ID0gQ2FzZXNQZXJNaWxsaW9uKSkgKwogIGdlb21fcG9pbnQoKSAKCgpgYGAKCkpPRVk6IEkgbGlrZSB0aGlzIGdyYXBoIGJlY2F1c2UgaXQgdmlzdWFsaXplcyBvbmUgb2YgdGhlIGZhY3RvcnMgdGhhdCBjb250cmlidXRlcyB0byB0aGUgZGV2ZWxvcG1lbnQgY2xhc3NpZmljYXRpb24uICBDYW4gd2UgYWRkIGluIGEgY29sb3IgY29kZSBvciBzb21ldGhpbmcgdGhhdCB3b3VsZCByZXByZXNlbnQgb25lIG1vcmUgZmFjdG9yLCBzdWNoIGFzIGRpZmZlcmVudCBjb3VudHJpZXMgb3IgcmVnaW9ucyBvZiB0aGUgd29ybGQuICBUaGUgbW9yZSB0aGF0IHdlIGNhbiBkbyB3aXRoIHRoaXMgb25lLCB0aGUgYmV0dGVyLgoKYGBge3J9CgpnZ3Bsb3QoZGF0YT1Hcm93dGgsIGFlcyh4PWhlYWx0aCwgeSA9IERlYXRoc1Blck1pbGxpb25NYXJjaEdyb3d0aCkpICsKICBnZW9tX3BvaW50KCkgCgoKYGBgCgpKT0VZOiAgSSBhbHNvIGxpa2VkIHRoaXMgZ3JhcGgsIGJ1dCB3ZSBuZWVkIHRvIHJlbW92ZSB0aGUgb3V0bGllciB0aGF0IHNjZXdzIGl0IGFuZCBtYWtlcyB0aGUgc2NhbGUgaGFyZCB0byBzZWUgYSByZWxhdGlvbnNoaXAuICBJdCBpcyB0aGUgb25lIHdpdGggYW4gWC12YWx1ZSBvZiBvdmVyIDEwLDAwMAoKYGBge3J9CmdncGxvdChkYXRhPUdyb3d0aCwgYWVzKHggPSBhaXJwb3J0cyx5ID0gQ2FzZXNQZXJNaWxsaW9uTWFyY2hHcm93dGgpKSArCiAgZ2VvbV9wb2ludCgpCgpgYGAKCkpPRVk6IEkgdGhpbmsgb25lIG1vcmUgZ3JhcGggaGVyZSB3aWxsIGJlIHN1ZmZpY2llbnQuIENhbiB3ZSBkbyBvbmUgdGhhdCBjb21wYXJlcyB0aGUgcmVsYXRpb25zaGlwIGJldHdlZW4gQ2FzZXMgcGVyIE1pbGxpb24sIERlYXRocyBwZXIgbWlsbGlvbiwgYW5kIHBvcHVsYXRpb24hISEgTWF5YmUgc29tZXRoaW5nIGxpa2UgYSBjbG9yb3BsZXRoIGdyYXBoIChJIHRoaW5rIHRoYXRzIHdoYXQgSSdtIHRoaW5raW5nIG9mKT8KYGBge3J9CgpgYGAKCg==